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NONPRINCIPAL-PLANE SCATTERING FROM FLAT PLATES 


SECOND -ORDER AND CORNER DIFFRACTIONS 


PART A 


NONPRINCIPAL-PLANE SCATTERING FROM FLAT PLATES 

SECOND-ORDER AND CORNER DIFFRACTIONS 

I . Introduction 

The modeling of a perfectly-conducting, rectangular plate for 
scattering in nonprincipal planes using the Method of Equivalent (MEC) 
currents was discussed in a previous report [1] . Two models using 
only first-order equivalent currents were presented. The first model 
used Geometrical Theory of Diffraction (GTD) equivalent currents 
[2] -[3], which are well behaved for monostatic scattering but contain 
singularities for bistatic scattering. A second model using Physical 
Optics (PO) [4] and Physical Theory of Diffraction (PTD) [5] 

equivalent currents was developed. These currents are well behaved 
for both monostatic and bistatic RCS predictions. The GTD and PO/PTD 
equivalent currents models give similar results and compare favorably 
with moment method (MM) and experimental results away from regions 
near and at grazing incidence. Near and at grazing incidence, 
higher-order scattering and corner diffraction mechanisms were thought 
to be significant factors in the total scattered field and a means of 
including these components was desired. 


* 


In this report , two new models of the plate for 
nonprincipal-plane scattering are explored. The first is a revised 
version of the PO/PTD model with second-order PTD equivalent currents 
[6} included to account for second-order interactions among the plate 
edges. The second model uses a heuristically derived corner 
diffraction coefficient [7], [8] to account for the corner scattering 
mechanism. The patterns obtained using the newer models are 
comparedto the data of the previously reported models , MM, and 
experimental results . 


ORIGINAL PAGE IS 
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II. Theory 




A. Second-Order PTD Equivalent Currents 

* * 

As with many versions of the MEC, the second-order PTD equivalent 
currents are formulated using the canonical perfectly conducting, 
infinite wedge geometry which is used to approximate other geometries. 
For the rectangular plate shown in Fig. 1, each edge is modeled as an 
infinite half-plane by setting the exterior wedge angle to 2rc . This 
is a valid approximation as long as the edges are electrically 
isolated; thus, the accuracy of this model increases as the electrical 
size of the plate increases. The general wedge geometry is shown in 
Fig. 2. This geometry is applicable to both first- and second-order 
equivalent currents. The directional vectors and angles are: 

A. 

s' - the unit vector in the direction of incidence. 

A 

s » the unit vector in the direction of observation. 

A 

t = the unit vector tangent to the edge of interest, 
directed so that it encircles the scatterer in a 
counterclockwise manner. 

A 

n = the unit vector normal to the edge of interest, lying 
on the illuminated face. 


fS 0 f = the angle between s' and the edge. 

A 

P 0 ** the angle between s and the edge. 

\jj r = the angle between the illuminated face and the 
edge-fixed plane of incidence. 

ip - the angle between the illuminated face and the 
edge-fixed plane of observation. 
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*■ > 


a = the skew angle of integration across the surface. 

In terms of the directional vectors, the directional angles are: 


/v 


r 


-1 * 

cos (s' • t) 


-1 

cos (s 


cos 


t) 


\f) = cos 


For monostatic scattering, \p - (//' , /3 0 



n - (3 ' , and s 


(la) 

(lb) 

(lc) 

(l d) 
For the 


flat plate, the unit vector in the direction of incidence for the 
geometry considered is: 

AAA A A 

s' = -a r ® -a x sinG'co s<p' - a y sinG'sin^' - a 2 cos©' (2) 

A different set of directional vectors and angles must be 
formulated for each edge. The vectors for each edge are: 


Edge 

1: 


= a y 

n l = “ a x 

(3a) 

Edge 

2: 

A 

t 2 

A 

“ “ a y 

n 2 = a x 

(3b) 

Edge 

3: 

A 

t 3 

A 

= " a x 

A A 

n 3 = _a y 

(3c) 

Edge 

4 : 

A 

A 

= a x 

A A 

” a y 

(3d) 


The resulting /3 Q ' functions for each edge are: 
Edge 1: cos/3 0 ' = -sin0'sin$' 


/ 2 2 
1 - sin 0' sin <p f 


Edge 2: cos/3 0 ^ = -cos/3 0 ' 

sin/3 0 ' = sin/3 0 ' 


(4a) 

(4b) 

(4c) 

(4d) 
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I 


Edge 3 : 

COS0 o ' 

«= sin0'cos$' 

(4e) 


sin0 o ' 

/ 2 2 
= /I - sin 0' cos <p ' 

(4f) 

Edge 4 : 

cos|3 0 ' 

- -cos0 o ' 

( 4g) 


sinp 0 ; 

“ sin0 o ' 

(4h) 

Using the /3 0 ' definitions, 

the ip' functions can be expressed as: 


Edge 1: 

COS 1//-L 

-sin0 f cos <P' 
sin0 o ' 

(5a) 


sin0 1 

COS0' 

sin f3 0 [ 

(5b) 

Edge 2 : 

COS0 2 

= -COSI//-L 

(5c) 


sin^ 

« sin^ 

<5d) 

Edge 3 : 

cos^ 3 

~sin0 f sin</>' 
sin/3 0 ' 

(5e) 


sirup 3 

cos0' 

sin0 o ' 

(5f ) 

Edge 4 : 

cosi p A 

= -COSlp 3 

(5g) 


sin \p A 

- sin i// 3 

(5h) 


For each edge the tangential components of the incident electric 
and magnetic fields are needed to determine the corresponding 
equivalent currents. Both soft and hard polarizations are considered. 
The incident fields are: 

Soft Polarization 

E l - J E„ (6a) 

H 1 - ; e (1/,) E„ (6b) 
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Hard Polarization 


i A -jk*r 

E = a 0 E 0 e 


<l/») 


(6c) 

H~ * -a j (1/7]) E^ e — (6d) 

To simplify the rectangular plate analysis, the incident fields 
are transformed to the rectangular coordinate system. The position 
vector, r, is : 


£ « a x x + a y y + a z z 

The propagation vector, k, for the incident field is: 

A A A A 

k *■ -k a r = -k(a x sin0'cos <p r + a y sin0'sin0 f + a 2 cos0' ) 


(7) 


( 8 ) 


With respect to the rectangular coordinate system, the incident fields 


are : 


Soft Polarization 

i jk (xsin0' cos<f>' +ysin0' sin<p' +zcos0' ) 

E — E 0 e 


H = (1/t?) E. e 


X [-a x sin$' +a y cos0' ] (9a) 

jk (xsin0' co s<p r +ysin0' sin <p' +zcos 0' ) 

a ^ 

X [a cos0' co s<f> r +a cos0' sin <f> f +a z sin0' ] (9b) 


Hard Polarization 


Ei - E 0 e 


jk (xsin0' cos #' +ysin0' sin#' +zcos0' ) 


^ = - (1/t)) E q e 


X [a x cos0' cos #' +a y cos 0' sin#' +a z sin0' ] 

jk (xsin0' cos#' +ysir>0' sin#' +zcos0' ) 


X [-a sin#' +a cos#' ] 


(9c) 


(9d) 
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The scattered fields for far-field observations are expressed in 


terms of the vector potentials : 

E r = 0 (10a) 

E 0 a -jw [A 0 + T)F^ ] (10b) 

E * 3 -5" IA * - ” F e 1 ll0cl 

The vector potentials are: 

/ - jkr 

A *= I(x' ,y',z') — dt' (11a) 

— 47T j — R 

C 

r -jkr 

F - ■— M(x' ,y',z') dV (lib) 

c 


For far-field analysis, the following simplifications can be made: 

/V 

R = r - r'cos# = r - r' • a r (for phase variation) (12a) 

R = r (for amplitude variation) (12b) 

where 

A /V 

r' = a x x + a y y (13) 

for the flat plate oriented as in Fig. 1. For phase variation: 

R = r - xsin0'cos$' - ysin0' sirup' (14) 

Finally, the vector potentials for far-field scattering from the 
plate become: 


A 


F 



c 


I (x, y, z) 


M (x, y, z) 


e 


jk (xsin0' cos(p ' +ysin0' sin0' ) 


dl 


jk {xsin0' cos<p f +ysin0' sin(f)' ) ^ 


(15a) 


(15b) 
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p 


The integrals are evaluated along the perimeter of the plate in a 


counterclockwise manner. On 

the surface of 

the plate. 

z=0 . Also, 

for each edge. 

the 

coordinates 

are : 



Edge 

1: 

x - a 

-b ^ y ^ b 


(16a) 

Edge 

2: 

x = -a 

-b ^ y ^ b 


(16b) 

Edge 

3: 

-a 3 x £ a 

y “ b 


(16c) 

Edge 

4 : 

-a ^ x ^ a 

y = -b 


(16d) 

To simplify 

the 

derivation, 

each vector 

potential 

integral is 


represented as a sum of four integrals , each corresponding to an edge 
of the plate. For the second-order currents, numerical integration is 
used in evaluating the integrals. DC and M must be determined 
separately for each edge taking into account the individual 
geometries. Again, each edge is viewed as the truncation of an 
infinite half-plane. 

The general geometry for the formulation of the second-order 
currents is shown in Fig. 3. The directional vectors and angles are 
for the edge of first-order diffraction defined above in the 

A 

description of Fig. 2. The diffracted ray travels along cr at a skew 
angle of fi 0 r . The axis of integration across the structure is again 

A 

cr, which effectively eliminates all singularities except the Ufimtsev 
singularity for forward observation at grazing incidence. The 
distance from the first point of diffraction, 0 X , to the second point, 
0 2 , is £. The tangent vectors to the edge of first diffraction is t x 
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Fig. 3. Geometry for second-order components. 


I 
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and the vector at the second edge is t 2 . These are oriented so that 
they encircle the scatterer in a counterclockwise manner. 

The second-order PTD equivalent currents components are 
defined as: 


m: 


T? 


2 2 2 

sin |3 2 


where 


*2 


2 

sin fS 2 


cr x t- 


[ 3 x K 2 ] 


t(t 2 X S) • K, ] 


-jk£ (O' 




-rf jk(T (O' * s) , 

j e d<r 


-rf 


(17a) 


(17b) 


(17c) 


The surface fringe current density, j , is expressed in terms of the 
exact fringe-current scattering solution to the wedge; therefore, K 2 
consists of a contour integral in the complex plane that is integrated 

a f 

along cr and evaluated only at the upper endpoint. is evaluated 
using the Method of Steepest Descent and the result for a half-plane 
is [5]: 


sf- 


AV~2 |<rx t 2 1 e 
2 

jksin /3 0 ' (fx+cos^' ) 


F [V L ] 


X 


V l-/x 


[n 1 (t 1 * H 0 ) cos (0' / 2 ) + H 0 )ficot^ 0 'cos (i/»' /2) 


+ (1/tj) (l-/i) csc/3 0 ' sin (^' /2) ] 


F [V L(l+cosi/> f ) 

V~2 


A 


[n l (t l* ilo) " (t l* H 0 > COt|3 0 'COS^' 
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0 


+ t x (t x * E*) (1/tj) csc0 o ' sini/i' ] 


(18) 


The expressions for the currents simplify to; 

f ~ f sinift, 

M 2 = ^ n 2 * £ 2 


f A f A f 

I 2 = t 2 * K 2 “ n 2 # ^2 COt/3 2 COSt// 2 


(19a) 


(19b) 


where K 2 is given in (18) for a half-plane. Other necessary 
quantities are: 

A A 

1 1# . t 2 * the unit vectors tangent to the edges of first- 
and second-order dif fractions , respectively, 

oriented so that they encircle the scatterer in a 
counterclockwise manner. 

A A 

n l' n 2 ** the unit vectors normal to the edges of first- and 
second-order dif f ractions, respectively, pointing 
inward and lying on the illuminated surface of the 
scatterer . 

A 

cr = the unit vector in the direction of integration skewed 
at an angle fi 0 ' with respect to the edge of the wedge 
so that it represents the grazing diffracted ray. 

or 


cr « n 1 sin/3 0 ' + t 1 cos/3 Q ' 


(20a) 


£ = the distance along cr from to 0 2 t the points of 
first- and second-order diffraction, respectively. 


or 


L = k i sin 2 /? 0 f 

F (x) = the modified Fresnel transition function. 
00 


(20b) 


F (x) = V j / Ti e 


jx 


/ 

J 


e ^ dt 


(20c) 
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1-2 


(20d) 


M 


cosy - cos 0 ' 


sin /3 ' 


sin 2 (y/2) 
sin 2 0 o ' 


y = the angle between <r and s . 

A A 

cosy « <r • s = sin/3 0 ' sin/^cos^ + cos/3 0 ' cosf^ 


(20e) 


P i = the polar angles of s in the coordinate system 
local to the edge of interest, the n i ,y,t i 
coordinate system. 


where 


A A 


cosp i = s • t i 

A 

A 

<20f) 

sin/^cosV'i = s 

A 

* n i 

A 

(2 Og ) 

sinp i sin^ i « s 

• y 

(20h) 


The directional angles (B 0 r , £ c , ip e , and i/j are defined in (4) and (5) . 

The second-order diffracted field is obtained by substituting 
and I 2 from (19) into the vector potential integrals of (11), 
integrating, and then substituting into (10) . Numerical integration 
must usually be used to find the integrals of (11) . The limits of 
integration on the integrals of (11) are found using ray tracing. The 
area that the first-order fields affect is bounded by the two extreme 
first-order diffracted rays. Fig. 4 illustrates this procedure. Edge 
AB illuminates the curve from A' to B' . Integration is along the 
boundary from A' to B' . Often illuminated regions overlap due to 
interactions from many edges. 

The total first- and second-order fields are found by adding the 
fields due to scattering by the PO currents, the PTD components, and 
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the second-order components. These fields are valid for all 
directions of illumination and observation, on and off the Keller cone 
of diffracted rays, except for the forward direction due to grazing 
incidence where an infinite singularity exists. This is the Ufimtsev 
singularity. 

In terms of the directional vectors, the directional angles are: 


COSp 2 ** s 

• t 



(21a) 

sinp 2 cos0 2 

A A 

= s • n 2 



(21b) 

sing 2 sin0 2 

A A 

= s • y 



(21c) 

For each edge of second-order 

diffraction. 

the directional vectors 

are : 





Edge 1 : 

A A 

t 2 l = a y 

A 

n 21 

A 

= " a x 

(22a) 

Edge 2 : 

A A 

t 22 = ~ a y 

A 

n 22 

A 

“ a x 

(22b) 

Edge 3: 

A /V 

t 23 = ~ a x 

A 

n 23 

A 

” a y 

(22c) 

Edge 4 : 

A A 

fc 24 = a x 

A 

n 24 

A 

a y 

( 226 .) 


Using the definitions of (22) in the equations of (21) along with the 
A sini p 2 cot/3 2 

definition of s' from (2), the factors sin/3 2 and co s\p 2 i n (19) 


simplify to the following for each edge: 

- , - sirups cost 

Edge 1: , ;L 21 = 

sm/3 21 


co s9' 


2 2 2 
sin Q' cos <p' +cos O' 


cot/3 21 cos^2i 


-sin 6' cos<p' sin0' 

2 2 2 ~ 
sin 0' cos (p' +cos 0' 


(23b) 
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Edge 2 : 


s . 2 „, 2 2 _ , 

“ sin 0' cos p'+cos 0' 


(23c) 


COtp 2 2 cos ^22 


2 

-sin 0' cos</>' sin0' 

2 2 2 

sin 0' cos <(>' +cos 0' 


(23d) 


Edge 3 : 


sin^ 2 3 

sin£ 2 3 


COS0' 

2 2 2 
sin 0'sin ^'+cos 0' 


(23e) 


C°tp23 co ®^23 


2 

sm 0'cos#'sin0' 

2 2 2 
sin 0' sin +cos 0' 


(23f) 


_ , . sinip ?A cosQ' 

Edge 4: , Zf * = 

sinfii* , 2 , 2 ,, 2 „ t 

sin 0'sm #'+cos 0' 


(23g) 


COt/3 24 COS024 “ 


sin 0' cos</>' sin<f>' 

2 2 2 
sin 0'sin <f>' +cos 0' 


(23h) 


The remaining factor in the equivalent currents equations that 
must be determined for each edge is K f 2 - The preceding vectors and 
functions refer to the edge of second-order diffraction only. The 
vectors and functions involved in the definition of involve the 

A A 

edge of first-order diffraction only except for the |<r X t 2 1 factor 
which involves directional vectors from the edges of both first- and 
second-order diffraction. For a half plane, K 2 is given in (18) . 

A /V 

The t 1 and n 1 vectors are defined for each edge in (3) . The 
incident fields, H 1 and E 1 are defined in (9). F (x) is the modified 
Fresnel transtion function of (20c) . The directional angles (3* and ip 


are the same as those introduced at the beginning of this section. 
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All the necessary functions of these angles derive easily from the 
relations of (4) and (5) . 

A 

The unit vector in the direction of integration, cr, for each edge 


Edge 1: 

A A A 

o'! = -a x sin£ 0 ' + a y cos0 o ' 

(24a) 

Edge 2 : 

A A A 

<r 2 ■ a x sin P o 2 " a y cos I3 q ' 2 

(24b) 

Edge 3 : 

a A A 

°3 “ _a x cos/3 0 ' - a y sin/3 0 ' 

(24c) 

Edge 4 : 

A A A 

°4 - a x cos/3 0 ' + a y sin0 o ' 

(24d) 


In general the term fi is defined as : 

2 

^ = cosy - cos <3 n ' 
sin 2 /3 0 ' 

where 


cosy = <r • s (26) 

For each edge, the y functions are: 

Edge 1: cosy^^ = -sin/3 0 ' sin0' cos<f>' +cos£ 0 ' sin0' sin0' (27a) 

Edge 2: cosy 2 = sin/3 0 ' sin0' cos<f>' -cos/3 o ^sin0' sirup' (27b) 

Edge 3: cosy 3 = -sin/3 0 ' sin0' cos#' -cos£ 0 ' sin0' sirup' (27c) 

Edge 4: cosy 4 = sin/3 0 ' sin0' cos^' +cos/3 0 ' sin0' sirup' (27d) 


Using these definitions and those of (4), the n' s for each edge 

are : 

cosy! - cos 2 /3 ' 

Edge 1: Mi = - (28a) 

sin 0 O ' 

cosy 2 - cos 2 /3 ' 

Edge 2: m 2 = (28b) 

sin / 3 0 ' 2 
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Edge 3 : fi 3 = 


Edge 4 : m 4 = 


cosy 3 - cos 0 O ' 

sin2 ^Q3 

2 

cosy 4 - cos |3 0 ' 

sin2 Po^ 


(28c) 


(28d) 


The coupling terms between the edge of first-order diffraction 

A A 

and the edge of second-order diffraction are the | <r X t 2 1 f the L, and 
the £ terms. For each pair of intersecting edges a different 
factor results due to these terms. Due to symmetry, one needs to 

o o , 

consider the range 0 :< <j>' ^ 90 only for the plate rotation angle. 

This eliminates having to consider interactions between some edges. 

The remaining interactions that one must consider are: 

(1) lst-order diffraction from edge 1 to edge 4. 

(2) lst-order diffraction from edge 1 to edge 2. 

(3) lst-order diffraction from edge 2 to edge 4. 

(4) lst-order diffraction from edge 2 to edge 1. 

(5) lst-order diffraction from edge 3 to edge 2. 

(6) lst-order diffraction from edge 3 to edge 4. 

(7) lst-order diffraction from edge 4 to edge 2. 

(8) lst-order diffraction from edge 4 to edge 3. 

A A 

The | cr X t 2 1 terms for each of these interactions reduce to: 


Edge 

1 

to 

Edge 

4 : 

Ki 

X 


= COS0 o ' 

(29a) 

Edge 

1 

to 

Edge 

2: 

A 

K 

X 

A 

“ sin P oi 

(29b) 

Edge 

2 

to 

Edge 

4: 

A 

K 

X 

fc a 4 1 

= COS0 o ' 

(29c) 
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c 


Edge 

2 

to 

Edge 

1 : 

1*2 

X 

t 2l i 

= sin|3 0 ' 

(29d) 

Edge 

3 

to 

Edge 

2 : 

A 

1*3 

X 

A 

fc 2 2 l 

= cos£ 0 ' 

(29e) 

Edge 

3 

to 

Edge 

4: 

A 

1*3 

X 

fc 2j 

= sin/V 3 

(29f ) 

Edge 

4 

to 

Edge 

2: 

A 

1*4 

X 

A 

t2 2 1 

= cos/3 0 ' 

(29g) 

Edge 

4 

to 

Edge 

3: 

A 

1*4 

X 

A 

fc 2 3 l 

= sin/3 0 ' 

(29h) 


The distance parameter, £, designates the distance from the point 
of first-order diffraction to the point of second-order diffraction 
measured along c r, the grazing diffracted ray. The £ parameters are 
constant functions of incidence angle only for opposite edge 
interactions but are functions of distance along the edge for adjacent 
edge interactions. Fig. 5 shows the geometry for determining the £ 
parameters for interactions between edges 1 and 2 and edges 1 and 4 . 
The geometries for the other interactions are similar. The limits of 
integration in the vector potential integrals of (15) vary according 
to the illumination of the edge of second-order diffraction. The 
extent of illumination is bounded by the grazing diffracted rays from 
the edge of first-order diffraction. Thus, the limits of integration 
are a function of the incidence angle. The distance parameters, and 
L, along with the limits of integration are given below for each pair 
of interacting edges. Recall that: 

I> = sin 2 /3 0 ' (30) 
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. > 


Edge _! to Edge _4 


- 


14 cos/3 0 r 


Ha ~ k ^i4 s i n /3 0 ^ 

Limits of integration: -a ^ x 3 a 


a - 2btan ) - x s a 


Edge !_ to Edge 2^ 


£ 


2a 


12 


c°s/3 0 ' 


L 12 — ^^12 s -*- n @o'j 

Limits of integration: 

-b £ y £ b - 


2a 


tan (jr-/3' ) 


Edge 2^ to Edge 4_ 


a + x. 


£■>* = 


24 sin/3 0 ' 

l 24 = k^ 14 sin 2 p o ; 


Limits of integration: 


-a £ x a 


for 0 O ' ^ 71 - tan (a/b) 
for 0 O ' 2: 7i - tan 1 (a/b) 


for 0 O ' ^ 71 - tan 1 (a/b) 


for 0 O ' £ tan 1 (a/b) 
-l 


-a £ x ^ -a + 2btan/3 0 ' for 0 O ' £ 7t - tan (a/b) 


Edge 2_ to Edge !_ 


£ 


2a 


' 21 sin/3 0 ' 


L 2i ~ ^^21 s ^ n Pqo 


Limits of integration: 
-b £ y ^ b - 


2a 


-1 


tan O') 


for 0 ' £ tan (a/b) 


] 
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(31a) 

(31b) 


(31c) 

Old) 


(31e) 
(31f ) 


(31g) 

(3lh) 


{ 3 1 i ) 


Edge 3_ to Edge 2 

b - y 2 
^ 32 ~ sin0 o ' 

l 32 = k^3 2 sin 2 ^ 0 ' 

Limits of integration: -b * y ^ b 


for p Q r £ tan (b/a) 


b - 2atan£ 0 ' 5 y < b for /3 0 ' £ tan 1 (b/a) 


Edge 3 to Edge £ 


sin 0o'. 


sin /3_ 


Limits of integration : 


a " Z — o-r - x < a 
tan /3 0 ' 


for |3 0 ' ^ tan (b/a) 


Edge _4 to Edge 2 

b + y 2 
^ 42 “ cos/3 0 ' 

L 42 - kje 42 sin 2 p c ' 

Limits of integration: -b ^ y 2 = b 


for /3 0 ' < 7C - tan (b/a) 


-b ^ y ^ -b + 2atan (7T-P 0 ' ) for £ 0 ' ^ tan 1 (b/a) 


Edge £ to Edge 3_ 


» = 2b 

43 COS0 o ' 

L 43 = k ^43 sin fio'i 


Limits of integration: 


3 " tan (ir-0 o ' ) X 3 


for 0 O ' s n - tan (b/a) 
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The second-order fields are determined by substituting the 
current components for each pair of interacting edges into the vector 
potential equations of (11) and using (10) . For opposite edges the 
integrals reduce to closed-form expressions. The adjacent-edge 
integrals must be evaluated numerically. Just as for the first-order 
dif f ractions, integration is with respect to either the x or y 
coordinate so that the integrals involved simplify to the following: 
Edge _1 to Edge 4 


I *= | e 
14 


j e F(V ) e 


jkxsin0' cos(f> f 


dx 


X 1 

X 


14 


f e ik- ^ 14 F (V l 14 (1+cos^ 1 ) ) e 


jkxsinG' cos <f)' 


dx 


Edge 1 to Edge 2 


L 12 


^ i2kysin0' sind)' 

j e dy 

*1 


exp 


= -2bsinc (2kbsin0' siruf >' ) - 


f 4ak ^ 1 

-j sin0'sin0' 

{ J tan (Tt-0 o ^) Y ) 

j2k sin0' sin$' 


(32a) 

(32b) 


(32c) 


Edge 2 to Edge _4 

T a $ "Dkifp * /— — 7- . jkxsinG' cosd' 

I 24 = J e ■ iq F(V L 24 (1-h 2 ) ) e dx (32d) 

x l 

4 m J e'^ 24 F(V L 2 7(1+cos 0 2 ) ) e jkxsin0 ' cos *' dx (32e) 

X 1 
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4 ak 


Edge 2 to Edge 1 


l 21 


^*7 j2kysin0' sin<p' 

J e dy 


yi 


exp 


2bsinc (2kbsin0' sin0' ) + 


b-; 


an Po 


— sin0' sin0' 


j2k sin 0' sirup* 


Edge 3 to Edge 2 


■Y 

yi 


I* = -f e ^32 F( /- L - 32 (l -M 3 ) ) e 3ky»i"0'-in^ dy 


T b - f 2 -jk£ 32 t* / / y 7 TT J — r V jkysin©' sin0' 

I = -j e F(v L 32 ( 1 +COS 0 3 ) ) e 

32 Yi 


Edge 3 to Edge 4 


l 34 


^ j2kxsin0' cos<p* 

J e dx 


exp 


( • 4 bk 

l 3 tan 0 o ; 


s in 0 ' cos<P* 


j2k sin0' cos<p' 


dy 


Edge 4 to Edge 2 


lj 2 = -f e' 3k ^2 F(V— ) e jkysine ' sin ^ 


-ik£. 


dy 


*1 


“jki? 42 . /— tt— ; 7 — — jkysinG' sin$' 

I 42 = -j e F(V L 42 (l+costfr 4 ) ) e J y ^ dy 


*1 


Edge _4 to Edge 3 


l 34 


- V 


j2kxsin0' sin <p* 


dx 


f . 4 kb . I 

exp 3 sin0'cos$' 

l tan ( 71 -/ 3 ') ^ J 


j2k sinG' cos<p' 


(32f ) 

(32g) 
( 32h) 

(32i) 

(32 j) 
(32k) 

(321) 
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B. Corner Diffraction 


Because the GTD and UTD diffraction coefficients are derived from 
the exact solution to scattering by an infinite wedge, the 
coefficients fail to account for the joining of two edges at a corner. 
For certain aspect angles, the scattering from the corners is 
significant. Pathak and Burnside developed a heuristic corner 
diffraction coefficient [7] -[8] based upon an appropriate, although 
non-rigorous, modification of an asymptotic evaluation of the 
radiation integral due to the equivalent edge current that would exist 
along the scattering edge if the corner were not present. One major 
flaw in the coefficient is that it is non-unique for certain 
backscatter angles near normal incidence [9], thus a 
rigorously-developed corner diffraction coefficient is desirable. 
However, Pathak and Burnside's coefficient is successful for many 
plate geometries and may be used with caution. 

The geometry for a corner in a planar surface is shown in Fig. 6. 
The total diffracted field from one corner is the sum of contributions 
from each of the edges comprising the corner. The general form of the 
corner diffracted field is: 


b‘(»)~e‘(o c ) .£ c /_ /5 


s+s c ) 


-jks 


(33) 


where 

E 1 ^) = the incident field at the corner. 

D c = the dyadic corner diffraction coefficient. 


- 28 - 




s = the distance from the corner to the receiver. 

s' = the distance from the source to the point of edge 
diffraction . 

s" - the distance from the point of edge diffraction to the 
receiver . 

s c = the distance from the source to the corner. 

~ i k s 

e « the phase factor. 

The dyadic corner diffraction coefficient, like the ordinary 
diffraction coefficient, is in terms of parallel and perpendicular 


components : 


~ A i A dc c A i A dc c 
D c= e h e h D h + e s e s D s 


(34) 


and Dg are the hard and soft corner diffraction coefficients. 


respectively, given by the following: 
- jir/4 


' s,h V 2nk Cs,h(QE) c°s/3 0C - cos/3 c 


y7in/3 0 sin/3 0 


F[kL c a(7r+/3 oc - 0 e >] 


(35) 


where 


-jTT/4 

c s ,h<Q E ) = x 

2v2nk sin/3 0 


r 

F[kLa (0 _ ) ] 

r La (0) /A 1 

cos (/3 12) 

< 

[kL c a (tt+|3 0C - 0 C )J 


T F [kLa 0 + ) ] 

[ La(0 + )/A 1 

> 

cos (0 + /2) 

[_kL c a (7t+0 oc - 0 C ) J 

J 


( 36 ) 
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F is the Fresnel transition function given by: 

00 

2 

F (x) = 2j v^x e jx J e dx (37) 

A 

+ 

are as: 

^ - 0 ± 0' (38) 

£ 0 is the Keller cone angle. The other angles and functions are: 


/3 C = the angle between the incident ray at the corner and 
the edge of interest. 

/3 oc = the angle between the diffracted ray at the corner 
and the extension of the edge of interest as shown in 
Fig. 3-5 . 


a (0) « 2cos (0/2) 


Q f O ** 2 

si " 


i . -ici 

c s c +s 


(39a) 

(39b) 

(39c) 


These fields simplify considerably for far-field, plane-wave 
scattering from the rectangular plate. 

The RCS of the flat plate in all planes can be determined using 
only the corner diffraction coefficient. Near and at normal 
incidence, this formulation fails due to the nonreciprocal nature of 
the corner diffraction coefficient. To alleviate this problem, near 
and at normal incidence, the GTD equivalent currents solution is used; 
and the corner diffraction results are used away from the problem 
area. The total scattered field from a corner consists of terms for 
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each of the adjoining edges; therefore, eight terms are needed for the 
rectangular plate. 

The general expression for the corner diffracted field due to one 
of the two adjoining edges is given in (33) . For far-field 
scattering, s = s c = s' = s” = oo ; therefore: ‘ 


E C (s) = -E 1 (Q c ) • D c — 


for far-field scattering. 

For far-field backscattering, L c ** oo so that: 
F[kL c a(7r+0 oc -0 c ) ] £ 1 
Also, /3 oc = n - p c and \p' - \p so that: 

V sing c sin/3 oc _ 
cos^ 0 - cos/3 2 


The diffraction coefficient simplifies to: 

- jTT/4 

<h = C Sfh (Q E ) (4 

v'Trck 

For far-field backscattering, L = L c = oo and C s h (Q E ) simplifies to: 


-e' j7r/4 f 1 

Cs,h(Q E ) = H ^ 

2V 2nk sin f3 Q [ ^ 2ticos <tt-/3 0 ) 


27rcos (n-/3 ) 


The total field scattered due to one edge adjoining at a corner 


simplifies to: 
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E i (Q c > 


j 

87tkcos/3 0 



X 



The top view of the plate geometry for the corner diffraction 
analysis is shown in Fig. 7. The incident fields are the same as in 
(9). The field must be determined at the corner of interest. The 
corners are designated as the following: 



Corner 

A: 

x = 

a 

y - 

b 






Corner 

B: 

X = 

-a 

y = 

b 






Corner 

C: 

X = 

-a 

y = 

-b 






Corner 

D: 

X = 

a 

y - 

-b 





The 

angles 0 are 

the 

same as 

those given in 

(5) . 

The 

0 angle 

used 

for 

the 

scattering component for a 

designated 

edge 

and 

corner 

is the 0 

associated with 

the 

edge . 

For 

example , 

the 

0! angle is 

used 

for 


scattering from corner A due to the presence of edge 1. 


The fi 0 angle is the angle between the -s' and the edge or: 

A 

cos/3 0 = -s ' • c 


(47) 


where 


c = the unit vector tangent to the edge of interest pointing 
outward from the corner of interest. 

The vectors and angular functions for each corner/edge combination 
are : 
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Corner A 


a /\ 



Edge 1 : 

C A1 

" a y 

cos /3 OA1 = • 

-sin0' sin<£' 

(48a) 


Edge 3 : 

A 

C A3 = 

A 

" a x 

c °sP 0A3 = ■ 

-sin0' cos <p r 

(48b) 

Corner B 








Edge 2 : 

A 

C B2 = 

A 

" a y 

c °s/ 3 0b2 = • 

-sin0' sin<p f 

(48c) 


Edge 3 : 

A 

C B3 “ 

A 

a x 

c°s| 3 0b3 - 

sin0' cos <p' 

(48d) 

Corner C 








Edge 2 : 

A 

C C2 “ 

A 

a y 

cos| 3 0c2 = 

sin0' sin$' 

(48e) 


Edge 4 : 

A 

C C4 = 

A 

a x 


sin0' cos <j>' 

( 4 8 f ) 

Corner D 








Edge 1 : 

A 

C D1 “ 

A 

a y 

cos ' 3 o d1 “ 

sin0' sin <p' 

( 48g) 


Edge 4 : 

A 

C D4 “ 

A 

a x 

cos| 3 0d4 = • 

-sin 0' cos 0' 

(48h) 

The final 

parameter 

that 

must be designated is s, the distance 

from 

the corner to the observation point. The 

geometry indicating 

these 

distances 

is in Fig 

. 8. 

For far-field scattering, the following 

approximations are used: 






Amplitude : 

s 

A “ 

1 S B = s c = S D = 

s 

(49a) 


Phase : 

S A = S 

- 

Si f ' / 2 (, 2 + b 

2 ) (cos</>' +sin0' ) 

(49b) 



s B = s 

- 

si f' J 2 (a 2 +b 

2 

) (sin$' -cos#' ) 

(49c) 



S c = S 

+ 

5i f' / 2(a 2 + b 

2 

) (cos<f>' +sin</>' ) 

(49d) 



s D = S 

- 

si f ' / 2 (a 2 +b 

2 

) (cos#' -sin#' ) 

(49f ) 


- 35 - 




I 

Edge 3 


Edge 


)f the plate for corner 


diffraction. 


The total field scattered by the plate consists of a term of the form 
of (46) for each edge joined at a corner. These eight terms are 
added to arrive at the total field. 
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III. Results 


Computations were made for a square plate with each side equal to 
5.73 X. Comparisons are made among the following: MM and 
experimental results , the PO equivalent currents model, the GTD 
equivalent currents model, the PO/PTD equivalent currents model, the 
PO/PTD/2nd-order equivalent currents model, and the corner diffraction 
coefficients model. 

o 

Fig. 9 shows soft polarization results for a 30 rotated plate. 
Even the PO equivalent currents, which account only for surface 
scattering, give good results near normal incidence. As the angle of 
incidence moves away from normal, there is a need for components to 
account for edge diffraction. The GTD and PO/PTD models greatly 
improve the results in the grazing regions, although there remains 
some disagreement which points to the necessity for higher-order and 
corner diffraction components. Fig. 10 shows results from the PO/PTD 
model with the second-order PTD coefficients added and also shows 
results from the corner diffraction coefficient model. Since the 
corner diffraction coefficients are inaccurate near and at normal 
incidence, the GTD equivalent currents solution is used in the region 

o 

±5 on either side of the normal direction and the corner diffraction 
solution is used elsewhere. The second-order currents do not improve 
the results for this case and even result in worse agreement than that 
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obtained using the first-order models. The corner diffraction model 
improves the results near grazing and agrees with the first-order 
models away from grazing, indicating that the corner diffraction 
mechanism is the more crucial scattering mechanism for this plate 
rotation and polarization. 

Fig. 11 shows the hard polarization results using the first-order 

o 

models for the same plate configuration rotated 30 . Near and at 

grazing there is a major discrepancy between the first-order, 
high-frequency models and the experimental and MM results. The PO/PTD 
models with the second-order components added and the corner 
diffraction coefficients model yield much better results although 
discrepancies still exist. These results are displayed in Fig. 12. 

o 

The soft polarization results for the same size plate rotated 45 
are in Fig. 13. One would expect that corner-diffraction would play a 
major role at this angle of rotation so that the large discrepancies 
of the first-order models near grazing incidence are not surprising. 
The addition of second-order equivalent currents yields excellent 
results in this case. Surprisingly, the corner diffraction model does 
not improve the results. These results, shown in Fig. 14, indicate 
that for this angle of rotation and polarization second-order 
components are the major contributing factors. 

Hard plarization results for the same rotation angle are shown in 
Figs, 15 and 16. Just as for the 30 rotated plate for hard 
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polarization, the first-order, high-frequency models do not produce 
good results near the grazing regions. The corner diffraction model 
greatly improves the results; however, the addition of second-order 
components does not result in an improved model. This indicates that 
the corner diffraction mechanism is dominant for this configuration. 

IV. Conclusions 

j 

The nonprincipal-plane scattering from a rectangular flat plate 
was considered. Comparisons among five high-frequency models, MM and 
experimental results were made. Near normal incidence all the models 

I 

agreed; however, near grazing incidence a need for higher-order 
and corner diffraction mechanisms was noted. In many instances the 

i 

second-order and corner-scattered fields formulated in this report 
improved the results. 

! 

I 


i 


- 48 - 


ORIQIWAt PA&L 
OF POOR QUALITY 


REFERENCES 


[1] C. A. Balanis and L. A. Polka, "Nonprincipal plane scattering of 
flat plates," Semiannual Report, Grant No. NAG-1-562, National 
Aeronautics and Space Administration, Langley Research Center, 
Hampton, VA, Jan. 31, 1989. 

[2] A. Michaeli, "Equivalent edge currents for arbitrary aspects of 
observation," IEEE Trans . Antennas Propagat., vol. AP-32, pp . 
252-258, Mar. 1984. 

[3] , "Correction to 'Equivalent edge currents for arbitrary 

aspects of observation'," IEEE Trans . Antennas Propagat., vol. 
AP-33, no. 2, p. 227, Feb. 1985. 

[4] , "Elimination of infinities in equivalent edge currents, 

part II: physical optics components," IEEE Trans . Antennas 
Propagat . , vol. AP-34, pp. 1034-1037, Aug. 1986. 

[ 5 ] , "Elimination of infinities in equivalent edge currents, 

part I: Fringe current components," IEEE Trans . Antennas 

Propagat. , vol. AP-34, pp. 912-918, July 1986. 

[6] , "Equivalent currents for second-order diffraction by the 

edges of perfectly conducting polygonal surfaces," IEEE Trans. 
Antennas Propagat . , vol. AP-35, pp. 183-190, Feb. 1987. 

[7] F. A. Sikta, W. D. Burnside, T. Chu, and L. Peters, Jr., 
"First-order equivalent current and corner diffraction 
scattering from flat-plate structures," IEEE Trans. Antennas 
Propagat ., vol. AP-31, pp. 584-589, July 1983. 

[8] P. H. Pathak, "Techniques for High Frequency Problems," in 

Antenna Handbook : Theory, Applications, and Design , Y. T. Lo, 

S. W. Lee, Eds. New York: Van Nostrand Reinhold, 1988. 

[9] A. Michaeli, "Comments on 'First-order equivalent current and 
corner diffraction scattering from flat plate structures,'" IEEE 
Trans. Antennas Propagat . , vol. AP-32, pp. 1011-1012, Sept. 

1984. 


- 49 - 


